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ABSTRACT 

Context. The Planck satellite, scheduled for launch in 2007, will produce a set of all sky maps in nine frequency bands spanning 
from 30 GHz to 857 GHz, with an unprecedented sensitivity and resolution. Planets, minor bodies and diffuse interplanetary 
dust will contribute to the (sub)mm sky emission observed by Planck, representing a source of foreground contamination to 
be removed before extracting the cosmological information. 

Aims. The aim of this paper is to assess the expected level of contamination in the survey of the forthcoming Planck mission. 
Metliods. Starting from existing far-infrared (far-IR) models of the Zodiacal Light Emission (ZLE), we present a new method 
to simulate the time-dependent level of contamination from ZLE at Planck frequencies. 

Results. We studied the possibility of Planck to detect and separate the ZLE contribution from the other astrophysical signals. 
Conclusions. We discuss the conditions in which Planck will be able to increase the existing information on the ZLE and IDP 
physical properties. 

This work is done in the framework of the Planck/Z/W activities. 
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1. Introduction 



The ESA Planck satellite ^ l)TaubeJl2003|) . scheduled for 
launch in 2007 ^, is a full-sky surveyor dedicated to cos- 
mic microwave background (CMB) and millimetric (mm) 
and sub-mm astronomy. It is a third generation microwave 
mission, after the NASA COBE and WMAP missions. The 
surveyor will observe the sky through a 1.5 m Gregorian 
aplanatic telescope carrying two instruments on the fo- 
cal surface operating at the frequency bands centred at 
30, 44, and 70 GHz (Low Frequency Instrument, LFI; 
iMandolesi erallllQQSll and 100, 143, 217, 3 53, 545, and 
857 G Hz (High Frequency Instrument, HFI; IPuget etlil] 
Il998|) . Planck will be injected in a Lissajous orbit around 
the Sun-Earth Lagrangian point L2 of the Sun-Earth sys- 
tem at a distance of 1507 683 Km ^ 0.01 AU from the 
Earth, from which it will observe the microwave sky for 
at least 15 months, necessary to complete two surveys of 
the whole sky with all the receivers. The LFI beams are 



Send offprint requests to: M. Maris 
^ http://www.rssd.esa.int/planck 

^ While this paper was at the end-stage of the editorial pro- 
cess, the launch has been posponed to 2008. This will not affect 
the results in this paper. 



located on the Planck telescope field of view in a ring 
with a radius of about 4° around the telescope line of sight 
(LOS) pointing at a scan angle a — 85° from the satellite 
spin axis. The HFI beams, located closer to the centre, 
may be also at few degrees from the LOS. In the simplest 
scanning strategy, the spin axis, chosen pointed in the op- 
posite direction from the Sun, will be kept parallel to the 
Sun-spacecraft direction, re-pointed by A^s = 2.5' once an 
hour (1° per day) in order to follow the revolution of the 
L2 Lagrangian point, and will spin at 1 RPM. Planck will 
scan the sky in nearly great circles approximately orthogo- 
nal to the ecliptic plane at a rate of 24 circles per day, each 
circle being scanned consecutively 60 times per hour. In 
this way Planck will produce at least two full sky maps 
for each frequency channel with an unprecedented resolu- 
tion (FWHM from ~ 33' to ~ 5') and sensitivity (in the 
range of ~ 10 — 49 mJy on a FWHM^ resolution element). 
Although the detailed scanning st rategy is currently un- 
der study l)DuDac fc TaubeJl2005(l . its general properties 
imply that only objects located at ^ 80° — 90° from the 
Sun will enter the large circles traced in the sky by the 
main beams. Solar System objects then will be observed 
nearly in quadrature with the Sun i.e. within few degrees 
from a solar elongation of ~ 85°. 
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Fig. 1. Relations between the hehocentric and cloud- 
centred frames used to describe the Solar System scan- 
ning geometry for a mission like Planck. The relations 
between Planck, Sun, Cloud Center and the observed por- 
tion of cloud are drawn, the L2 point is not drawn to sim- 
plify the graph. The connection between the other two 
reference frames is shown in Fig. [3 An example of scan 
circle and the related spin axis, S, are also drawn. The 
graph is not in scale with real distances. (Colour figure on 
the electronic version.) 



Why would a cosmological mission like Planck con- 
sider the Solar System in its scientific program? How 
should Solar System studies take advantage from a mis- 
sion like Planck? A valuable contribution has been 
provided by the CO BE mission to solar system studies 
^Kelsall et al. 11X99^ . The scanning strategy of Planck 
assures that all the Solar System components located out- 
side the Earth Orbit, except Mars, will enter at least once 
into the field of view of the surveyor during the mission. 



etR:lBiirie-ana et al.'"200l' .'Pasre et al.''2003'. 
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comets) and the thermal emission of the diffuse in- 



terplanetary dust will represent foregrounds which have 
to be detected and properly removed in order to avoid 
the introd uction of systematic errors in the cos mological 
measures (jMaris. et al.ll2003t iMaris et al]l2004a|) . 

The Zodiacal Light Emission (ZLE) due to thermal 
emission from the Interplanetary Dust Particles (IDPs) 
is the far-IR counterpart of the familiar Zodiacal Light 
due to scattering of the solar light by IDPs. Most of the 
proper ties of the ZLE below 300 /um have been studied by 
IRAS fWheelock et al.'l994'), COBE fKelsall et al. "l998|) 

and ISO (R each et al..Ll99fi; .Reach et al 200.1') . Peaking 

at A « 10 /im, the ZLE is one of the major contributors to 
the sky background in the far-IR domain at low ecliptic 
latitudes. A first dete ction in the 300 um — 1 000 iim band 
has been assessed by iFixsen fc Dwekl l|2002() using yearly 
averaged COBE/FI RAS data. Even a q uick look at the 
data r eported in .Kelsall et al. I l)l998l) and lFixsen fc Dwekl 
( 2002() allows us to predict a contribution of ZLE in the 



857 GHz channel of Planck of « 0.6 MJy/sr. It is evident 
that at the Planck lowest frequencies its contribution 
is much weaker than the Galaxy emission. On the other 
hand, at the Planck intermediate and high frequencies 
the ZLE is significantly weaker than the Galactic emission 
only at low Galactic latitudes while it is comparable to it 
outside the Galactic plane (for example, near the poles 
the Galactic emission is ~ 1 MJy/sr at 857 GHz). The 
expected ZLE contribution is larger than the instrumen- 
tal noise at the Planck highest frequencies. Therefore, 
a careful analysis of the ZLE in the Planck data is re- 
quired. Since the ZLE varies over angular scales ^ 10° 
it can be properly studied by working at a resolution 
of ~ 1° — 2°. At this scale the expected la sensitivity 
per FWHM^ at 857 GHz at the e nd of the mission is 
w 2 X 10"^ MJy/sr ijLamarre et al.1.200 3). The extrapo- 
lated background from the Galaxy, representing the main 
large scale background component at this frequency, is 
^ 1 MJy/sr. Because of the different tilt on the ecliptic of 

the Galactic plane and of the IDP cloud, for most of the 
scan circles observed by Planck, the sky position of the 
maximum of the ZLE will fall close to that of the minimum 
of the Galactic emission. In this case the ZLE is extrapo- 
lated to be about half of the Galactic emission. Of course 
one may wonder whether ground-based or balloon born 
CMB experiments may have been able to detect such a 
contribution. Looking at some of th e most recent ball oon- 
borne experiments, ARCHEOPS (Benoit et al.ll2003l) has 
constructed maps by bandpassing the data between 0.3 
and 45 Hz, correspond ing to about 30° and 15' scales, re- 
spectively. MAXIMA ijLee et al.ll200ll) covered the mul- 
tipole range 36 < i < 1235 or angular scales smaller 
that 6°. The CMB po wer spectrum of BOOMERANG 
llNetterfield et alJl2002^ covers multipoles from ^ ~ 75 to 
£ ~ 1025, equivalent to angular scales 10' — 2.4° From 
ground based experiments, DASI ijPrvke et al.l200^ mea- 
sured the power spectrum for 100 < £ < 900 or angular 
scales less than 1.8°. Among these experiments it seems 
that only ARCHEOPS is able to detect the large scale 
brightness variation connected with the ZLE. We then ex- 
pect that the ZLE will be observable as an excess of signal 
superimposed on the Galaxy or it may be considered as a 
source of systematics in studying the large scale Galactic 
emission. 

Our strategy, in line with past studies, exploits the ex- 
isting far-IR observations, included in models, to derive 
the spatial distribution of the ZLE and to extrapolate its 
Spectral Energy Distribution (SED) at Planck frequen- 
cies. We take as a refe rence the yearly averaged values of 
iFixsen &: Dwekl l|2002j) . Our starting point to model the 
spatia l distribution of the ZLE, the work of iKelsall et al. I 
(1998^ for the ZLE based on the COBE data (hereafter in- 
dicated as the C Q BE-model), has many s imilarities with 
the IRAS model bv lWheelock et alJ l)l994|) . It describes in 
detail the emissivity of the IDP cloud, assumed to extend 
up to ~ 5.2 AU from the Sun, for wavelengths up to about 
300 /Ltm. According to the COBE-model four components 
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Fig. 2. Relations between the cloud-centred ecliptic ref- 
erence frame, the cloud symmetry reference frame and 
the spacecraft centred corotating reference frame. Top 
panel: angles between the cloud-centred Ecliptical refer- 
ence frame (full yellow disc (light-gray in the bw version) 
- black arrows), and the cloud cylindrical symmetry refer- 
ence frame (dashed blue disc - red arrows). Bottom panel: 
angles between the cloud cylindrical symmetry reference 
frame, the cloud-centred corotating cylindrical reference 
frame, the spacecraft - centred corotating reference frame. 
The displacement of Planck (in this case above the sym- 
metry plane), ^Rp, respect to the constant reference po- 
sition, Rp, assumed to be in the symmetry plane of the 
corotating reference frame is also displayed. The graph 
is not in scale with real distances. (Colour figure on the 
electronic version.) 

contribute to the ZLE: the dominating smooth compo- 
nent, the Earth orbit locked ring of dust (or circumsolar 
ring), the trailing blob, and three bands of dust. 

In this work only the standard IDP component has 
been considered. The analysis of the plausible, but not yet 
determined, cont ribution from the Kuiper Belt dust grains 
l|Landgraf et al.ll2002.) will be the subject of another work. 

With respect to other foregrounds usually considered 
in CMB studies, the ZLE (as the other Solar System ob- 
jects) is peculiar, depending for its surface brightness not 
only on the pointing direction but also on the instanta- 



neous position of the observer within the Solar System. 
Galactic and extragalactic foregrounds are generated by 
sources located so far from the observer that parallactic 
effects due to the motion of the observer within the Solar 
System are negligible compared to the instrumental res- 
olution of CMB observatories. On the contrary, since the 
observer is located within the Solar System, the orbital 
motion about the Sun leads to changes in the ZLE bright- 
ness distribution as a function of the pointing direction. 
This underlines the relevance of studying the ZLE not only 
on maps but also on time ordered data streams (TODs). 
Moreover, the accurate simulation of observations for a 
satellite mission like Planck cannot be based on maps 
since the details of the orbit will have to be considered in 
addition to the usual scanning law. 

The main aim of this work is to contribute to the fol- 
lowing subjects: i) to define a representation method for 
the ZLE suitable for map based CMB mission simulators, 
with particular relevance for the Planck mission; ii) to 
determine suitable approximations (like scaling frequency 
laws) for the simulation of this component in the frame- 
work of the Planck simulation pipeline; iii) to determine 
to what extent the ZLE will impact the Planck survey; 

iv) to determine if it will be possible to separate the contri- 
bution of the ZLE from the data produced by the Planck 
mission in a self-consistent manner (i.e. reducing as much 
as possible the introduction of priors based on results from 
other missions in the Planck data processing pipeline); 

v) to explore the possibility of Planck to produce use- 
ful scientific results about the ZLE at frequencies barely 
explored in the past. 

The paper is organised as follows. In Sect.|21we briefly 
review the present knowledge about the ZLE. Sect. 01 de- 
scribes the fram ework of our s imula tions, mainly based 
on the model of iKelsall et al. I ((l998^ . and the details of 
our numerical code discussing its main assumptions in the 
light of recent theoretical results. In Sect. I3.ll we present 
a series expansion of the ZLE spatial dependence that can 
be useful for many simulations and data analysis applica- 
tions. In Sect.^we describe the main results of these sim- 
ulations (mainly in form of TODs and maps) and compare 
the ZLE contribution to those expected from the Galactic 
emission. Sect. |5l is devoted to the separation of ZLE in 
the Planck data. Particular care is given to the analy- 
sis of the systematic effects in the differential approach 
for ZLE separation. Our main results and conclusions are 
summarised in Sect.El 

2. Physical and geometrical properties of ZLE 

To assess the expected errors in predicting the ZLE sur- 
face brightness at Planck frequencies from a model based 
on far-IR and IR data we review some theoretical con- 
cepts needed to link the ZLE model, and in particular 
the COBE-model, to the optical properties and the size 
distribution and other physical properties of the IDPs. In 
sub- mm and mm bands the dominant emission mechanism 
from IDPs is thermal emission of IR radiation driven by 
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solar heating. The most general expression for the bright- 
ness averaged over the bandwidth and the beam, detected 
in a radiometric channel of frequency / on-board a space- 
born experiment produced by a given population of IDP 
grains, representing a component c of the IDPs cloud is 



/y,,(P,Rp) = 



/ dvdP'Wfiv) Gf(P',iy) 
Jn 



(1) 



Jo Jo 

B,(r,(r(s))7ra2Q,b.c(i^, r(s)) ^^^^^^ 



here a is the grain size, i/ the frequency, P is an ob- 
serving direction, s the distance from the observer along 
P, Rp the position of the observer respect to the Sun, 
r(s) = sP -I- Rp the position within the Solar System 
respect to the Sun along the line of sight, W{i') is the 
instrumental frequency response, G^{P) the beam re- 
sponse, Tc(r) is the thermodynamical temperature for 
grains in the population c, B,y the blackbody emissiv- 
ity and Qnhsd^i o,,^) is the absorption coefRcient for the 
population c of grains, dNc{a,T:)/da the size distribution 
for grains. A rigorous calculation of this integral is diffi- 
cult since it requires the knowledge of many poorly con- 
strained quantities. Alternatively effective models, such 
as the COBE model, are used in place of Eq. We dis- 
cuss here the relation between the two approaches. Let us 
assume that: i) the size distribution and the grain optical 
properties do not depend on r, QahsA^, a, r) ss Qahsd'^, a), 
dNc{a,r)/da « dnc{a)/da ■ Nc{r) with Nc{r) the spatial 
distribution of grains; ii) the beam is symmetrical and it 
does not depend on v within the frequency bandwidth of 
each channel /; iii) the beam is small compared to the 
typical angular scales over which the ZLE varies. With 
these assumptions Eq. simplifies to 



//,c(P,Rp) = KfEf,cno,c 

./+°°d.i?/(r,(r(.))iV,(r(.)), 



(2) 



where Kf is a color correction which takes into account 
the averaging over the frequency weighted by the instru- 
mental response (see Sect. l2.2|l .. no,c is the optical density 
for the given dust component, Ef,c is an emissivity cor- 
rection related to the size distribution and composition 
of grains. For multifrequency observations, Ef is usually 
normalised to a refe rence frequency, /q. For example in 
iKelsaJ] et al.l ^98) Ef is normalised to the value it has 
in the COBE/DIRBE 25 /xm (/o = 12 000 GHz) chan- 
nel. Note that such kind of normalisation Ef implies a 
corresponding renormalisation of the optical density no,c 
in Eq. In principle it would be possible to compute 
Qabs(f^7 a) and then Ef from Mie theory assuming appro- 
priate grain shapes and composition and usin g commonly 
available software jBohren fc Huffmanljl99^. as done re- 
cently in the studv bv iReach er'aLT(f2003|) of the ZLE 
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Table 1. 



-^/^Smooth 



estimated according to Eq. 



emissivity at A < 100 fiia. However, the composition and 
shape of IDP grains in the size range relevant at Planck 
frequencies are only poorly known. We prefer to deter- 
mine Ef at Planck frequencies comparing theoretical es- 
timates with the existing COBE/FIRAS data, as detailed 
in the next section. 



2.1. Estimating Ef at Planck frequencies 

We estimated Ef in the relevant range of frequencies com- 
paring the expected ZLE yearly averaged from the COBE 
model with the existing surface bri ghtness measures at 
Planck frequencies from lFixsen fc Dwek. (,2002.) based on 
COBE/DIRBE and COBE/FIRAS data. However, given 
the uncertainties in the interpretation of these data, other 
extrapola tion methods are possib le. According to the main 
result of iFixsen fc D^ l)2002|) . the SED of the ZLE 
Smooth component is approximately similar to a black- 
body with T « 240 K scaled by an emissivity factor nearly 
constant for A < 160 /im and scaling as A~^ at longer 
wavelengths: 



(^.=/)y£r « 3 X 10-^ X B,{T - 240 K) 



1875 GHz; 



> 1875 GHz, 
' , 1/ < 1875 GHz 



(3) 



As evident from this analysis of COBE/FIRAS data 
llFixsen fc Dwekl |2002|) . the accuracy of the recovered 
ZLE spectrum is good at frequencies higher than about 
800 GHz while the error bars significantly increase at lower 
frequencies, where most of the Planck frequency chan- 
nels are located. This leaves space for further improvement 
in the field with Planck. After estimating the yearly av- 
erage surface brightness from the COBE model, (Z/}year, 
assuming Ef = I we obtain 



E 



FD 

/, Smooth 



/ 7- \FD 
\^J'=//yoar 



(4) 



year 



The resulting estimated Ef, normalized to 1 at / = 1.2 x 
10"* GHz, are reported in Table^with the expected yearly 
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average surface brightnesses. E^^^^^^^^^ is approximated in 
the range of frequencies of interest by 



Smooth, 857 GHz 



MJy/sterad 



E 



FD 

/, Smooth 



126.63 



/ 



1.2 X lO^GHz 



(5) 



2.2. The colour correction 



The colour correction Kj in Eq. Q takes into account 
the effect of the frequency instrumental response within 
the bandwidth. At Planck frequencies, most of the ZLE 
is due to IDPs with temperatures exceeding « 200 K; 
within the bandwidth of each Planck frequency chan- 
nel, their emission integrated along the line of sight can 
be approximated by a power law with spectral index 7/ 
and normalisation J^f{v = f) ~ Ff (from here we will 
omit the pointing dependence to simplify the notation). 
In addition, within a frequency band Tf(y) = €f{v)Zf{v), 
where e/(i/) = (yjj'yi-'Ef is the spectral emissivity cor- 
rection and Zf{v) = {v/ f)^f-'Zf accounts for the spatial 
distribution expected from grains emitting as blackbod- 
ies and can be separated into a pure spatial dependence 
Zf{v) and a pure frequency scaling. Here normalisations 
are defined to match the values of Ef and Zf ior v — f 
and the overall 7/ turns to be 7/ = ^f^^ +Jf,z- Therefore, 
by imposing 



FfEfKf^^f Wf{v)e{v)Tf{v)dv 
Jew 



(6) 



and Taylor expanding in the log - log space about log/ 
the brightness and the emissivity correction, we obtain 



1 



^/ Jbw 



dv 



(7) 



together with Ff = Tf{f) and Ef = e(/). The validity of 
the approximation represented by Eq. (0) is then verified. 

We assess the value of Kf needed to compensate the 
fact that in Eq. Q the integration over the bandwidth in 
Eq. has been neglected. In addition we want to assess 
the level of uncertainty in the Kf correction induced by 
the uncertainty in the Ef prediction. 

We restrict ourselves to the illustrative case of a sim- 
ple top- hat window Wf{y) with relative bandwidth rf — 
Af/f. For Planck r/ = 0.2 for / < 70 GHz and 
r/ = 0.25 for / > 100 GHz. We tested that for differ- 
ent reasonable shapes of Wf (ly) the results do not change 
significantly. Under these conditions, for 7/ 7^ 1 



Kf 



r/(l -1-7/) 



H7/ 



T-y-j 1+7/ 



(8) 



where Wqj is the normalisation for the spectral window. 

For numerical estimates, assuming Wqj = 1 and 
7/,z ~ 7/,e ~ 2, we obtain A'/<7o ghz = 1.020 and 
-f''/>ioo GHz = 1.031. For the case of a frequency- 
independent Ef, 7/^(: = and we find i^/<7o GHz = 1.003 
and iC/>iooGHz — 1.005. Leaving 7/ to vary within ±1 




-100 -50 50 100 

Relative Ecliptical Longitude [deg] 



Fig. 3. Simulated contour plot of Zf{P) [MJy/sr] for a 
fixed location within the Solar System. Contours are cal- 
culated for Zf = 0.205, 0.22, 0.23, 0.24 MJy/sr and from 
0.25 MJy/sr up to 3.2 MJy/sr in steps of 0.05 MJy/sr. 
Surface brightnesses exceeding 3.2 MJy/sr occurring in 
the central empty region have been neglected. The blue- 
dotted line represents the path of a scan circle of Planck 
for an horn located at the centre of the field of view while 
the gray band represents the Planck field of view. The 
vertical axis is the latitude over the ecliptic. The horizon- 
tal axis the longitude relative to the solar direction (oppo- 
site to orientation of the scan axis) . Note that in this case 
the symmetry plane of the IDP cloud is below the ecliptic 
plane. (Colour figure on the electronic version.) 

unit Kf changes by only 2%. Thus it will be possible to 
avoid applying this small correction in the numerical esti- 
mates presente d in the remaining pa rt of this paper. 

Note that in iKelsall et al. I l(l998l) the colour correction 
is defined as a correction for the instrumental response 
when a blackbody is observed. For this reason the colour 
correction is parametrised as a function of the blackbody 
temperature T which is a function of the position along the 
line of sight, so that the argument of the pointing direction 
integral in Eq. |5J) would have to be scaled by the spatial 
dependent Kf{T). However, at our frequencies the bulk of 
the blackbody emissivity comes from grains emitting not 
too far from the Rayleigh- Jeans limit, i.e. with a frequency 
power law scaling (within our rather limited bandwidth) 
with a power law index largely independent of T. In this 
case the two definitions for Kf are equivalent. 

3. The model and the numerical code 

The surface brightness calculated for a given frequency 
band is 



-F'zL_B,/(P,R0,R-p) — Z]c-^c,/(P,R0,R-p) 

= J2c EcjZcj(P, Rq, Rp) . 



(9) 



Here R© denotes the position of the Sun within the Solar 
System and the index c denotes the specific component 
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either the dominating smooth component, the Earth or- 
bit locked ring of dust, the trailing blob, or one of the 
three bands of dust. According to the discussion in Sect.|21 
in Eq. @ we separate the calculation of the spatial dis- 
tribution of the ZLE assumed to be a blackbody from 
the more uncertain emissivity correction. The detailed 
COBE MODEL has been already describ ed in the litera- 
ture l|Ke]sa,]l et aL IPixsen fc n^^l2nn^ and only 
few details need to be reviewed. The relevant geometry 
of ZLE observations is shown in Fig. ^ the Sun (or the 
barycentre of the Solar System Bss), the position of the L2 
point, the spacecraft (in this case Planck) position Rp, 
the centre of the distribution of dust related to the com- 
ponent c of interest, Cc, and the related vectors drawn 
between these points. The most important among them 
are: the position of a point at distance s from the space- 
craft along the pointing direction P: R(s) = Rp + sP; its 
position with respect to the Sun, Rg = R — Rq, and to 
the centre of the cloud defined by Rq^c- R'c — ^ ~ R-o,c- 
In addition we define R'^ — |R^|. X'^, Y^, Z'^ denote the 
Cartesian components of R^ and X, Y, Z those of R. 
The model specifies the 3D dust density distribution of 
each component, factorised into a radial dependence and 
a vertical dependence, as the heliocentric dependence of 
the mean dust temperature. In this way, for each com- 
ponent c it is possible to define a proper reference frame 
about which the cloud c has cylindrical symmetry. Its ori- 
gin coincides with the centre of the cloud and its midplane 
with the equatorial plane of the cloud; it may be rotated, 
tilted and shifted with respect to the ecliptic plane. 

We have implemented the COBE MODEL in 
a FORTRAN-90/95 program called FS_ZDD (FLIGHT 
Simulator - Zodiacal Light Emission) embedded in 
a supporting OCTAVE pipehne fjMarisil2001lV The code was 
originally designed as a module of the Planck Flight 
Simulator but can be used for any other experiment. Since 
the Ef scalings are largely uncertain it is left to the user to 
apply the proper one to the output of these programs. In 
this work we use the code to study the time dependence in 
the signa l acquired by Planck for the nominal scanning 
strategy llDupac &: Taubeil2005l) and a recently simulated 
spacecraft orbit jHechle^2002l) and to predict the ZLE in- 
duced perturbation in Planck data. We choose to express 
our results in terms of brightness (MJy/sr), as for IRAS 
and COBE data. 

3.1. Series expansion of the ZLE spatial dependence 

Simulators for CMB missions and the related data- 
reduction pipelines, such as those reahzed for Planck, are 
largely based on maps. A map allows a good representa- 
tion of the sky brightness as a function of the pointing di- 
rection. This procedure neglects the time-dependent infor- 
mation on the Planck position within the Solar System, 
leading to a loss of information when Solar System com- 
ponents are considered. The ground segment of a mission 
like Planck would be able to handle and analyse TODs as 
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Fig. 4. Simulated data stream of surface brightnesses 
(MJy/sr) measured at 857 GHz for the ZLE - smooth 
component (red), the Galaxy (green) and the sum of the 
two (blue). The ordinate is the phase of the scan circle, 
assumed to be zero for the pointing direction nearest to 
the North ecliptic Pole. Two subsequent scan circles are 
displayed, the phase of the second being augmented of 
360°. The lowest, full, cyan line (peaks at 0.1 MJy/sr) 
represents the contribution of secondary components of 
the ZLE. The red-dashed line represents the full contribu- 
tion of ZLE with both smooth component and secondary 
components. The cyan band represents the uncertainty 
in the prediction for the Smooth component due to the 
uncertainty in the extrapolated Ef value in the case of 
TOc ^ 10~^ gr. Both the dashed line and the cyan band are 
reported after adding the Galactic contribution as a green 
dashed line and a grey band above the Galaxy. (Colour 
figure on the electronic version.) 

well as maps obtained from them (|P asian fc Svgnet"2002'; 
Ivan Leeuwen et al.l200ilChallinor et al..2002tl . A TOD of 
ZLE would allow an exact representation of any seasonal 
dependence. However, TODs are large and their realiza- 
tion requires an effective scanning strategy and satellite 
orbit, possibly accommodated during the mission, so that 
the exchange of simulated data in the form of TODs is 
not practical for the data analysis of a multichannel, high- 
resolution mission like Planck. We implement a method 
able to 1. properly represent seasonal effects in a large 
set of mission configurations; 2. possibly be applied to 
other missions; 3. exploit the (cylindrical) symmetries in 
the components of the IDP cloud. 

We propose to generalise the concept of pixelized map. 
A pixelized map is usually defined as the values assumed 
by a given observable on a set of pixels ordered according 
to the adopted pixelization scheme. In this contest we can 
introduce a generalisation of this concept by defining a 
pixelized "map" as a list of values assumed by a given ob- 
servable on a set of pixels which are also functions of the 
positions of the Sun (R©) and of the Spacecraft (Planck 
in our case) (Rp) within the Solar System. In a reference 
frame (r.f.) in which the displacements of the Sun and 
spacecraft positions are just small fractions of their av- 
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eraged positions, the pixehzed "map" can be replaced by 
the series expansion of the observable about the average 
positions of the Sun and the Spacecraft. 

Then, denoting with Rp, the reference positions 
of Spacecraft and Sun about which the series expansion is 
performed, we have 



57 GHz 



Rp — Rp 



^Rf 



Rm — R(r 



SKr. 



(10) 



In the numerical computations for Planck we exploit val- 
ues of the displacements along each direction up to about 
±0.03 AU for SKq and up to about ±0.07 AU for (5Rp, 
about a factor of two wider than any reasonable displace- 
ment for a mission that will reside near h^, . In addition, 
we ad opt in this work the HEALPix scheme I Gorskv et al.l 
120051) with "ring" ordering, widely used in the CMB com- 
munity. 

We have tested that in the case of the ZLE it is 
preferable to expand not the brightness spatial distri- 
bution but its logarithm in power series. Then, denot- 
ing by Z/_c,p(R-0, R-p) the brightness integral for a given 
frequency channel /, component c and pixel index p 
(connected to the pointing direction P by the mapping 
scheme) as a function of R© and Rp we adopt the follow- 
ing decomposition 



^/,c,p(R-0, Rp) — ^/,c,p^s,/,c,p('5R0)^p,/,c,p(^R-p) 
•-Fps,/,c,p(^R-0, <5Rp), 



(11) 



where Zf^c,p = 2'/^c,p(R.0, Rp) and Fsj^cpiSRo), 
-Pp,/,c,p(<5Rp), ^PS,/.c.p('^R0, <5Rp) are exponential func- 
tions of polynomials of (5Rq and (5Rp . In order to achieve 
an accuracy better than ~ 1% for more than 96% of the 
pixels, we verified that ^ it is sufficient to consider the 
following terms of the series expansion: 



logFs,p(<5RQ) 
logi^sP,p('5R0<5Rp) 



as,p,iSRQ^i (12) 
+bs,p,ij5RQ,iSRQj , 

,SRp^, (13) 
+bp^p,ijSRp^iSRpj 

+ Cp,p,ijkSRp,iSRpjdRp^k , 

bsp,p,ijSRQ,iSRpj , (14) 



where SRQ^i and SRp^i, i = X,Y, Z are the Cartesian com- 
ponents of SRq and (5Rp; the a, b, and c terms denote the 
first, second and third order coefficients, respectively, and 
repeated indices i, j, k = X , Y , Z are summed. The num- 
ber of independent coefficients is 3 for ag.p^i and ap,p^i, 6 
for bs,p,ij and 6s,p,ij , 9 for cp.p^iju and 6sp,p,y - The other 
coefficients are simply obtained by index permutations. 
Then, for each pointing direction a total of 38 independent 



^ The remaining pixels will have an accuracy worse than 1% 
(the worst accuracy recorded being 2.5%) in the case of a dis- 
placement of the spacecraft and of the Sun at the limits of the 
region for which the expansion is calculated. 
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Fig. 5. Relative contribution of ZLE, Galactic (dust) 
emission and noise as a function of the pointing eclipti- 
cal longitude at three Planck frequency bands: 857 GHz 
(top), 545 GHz (middle), 343 GHz (bottom). For each 
pointing ecliptical longitude, the plots gaves the ratio of 
the ZLE over the Galactic emission EfZf/Gf averaged 
over the given scan circle (white - full line), the ±3(T range 
(blue (black in the bw version) band) , the peak ratio (yel- 
low (gray) - full line) and the averaged ratio between the 
instrumental noise and the Galactic emission. Data are 
calculated for patches of 1° in radius, noise is for a 14 
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Fig. 6. Derivatives of ^smooth, 857 GHz along (blue full line 
- (5||) and across (red dashed line - 5±) the scan circle of 
Fig. 121 For graphical reasons derivatives in the plot are 
scaled by a factor 10'^. According to the convention of 
Fig. ^ the phase represents the pointing position along 
the scan circle, with for the pointing direction nearest 
to the North ecliptical Pole. Note that the 6± direction is 
always normal to and oriented toward the local solar 
direction so that din Zf /dS± is always positive. (Colour 
figure on the electronic version.) 

components has to be computed. These coefficients can be 
determined by solving by least squares a sufficiently large 
set of independent equations obtained by exploiting differ- 
ent combinations of displacements (5Rp , SHq (in the case 
of Planck about 100). Metrics are then applied to assess 
the quality of the data generated with this series expansion 
compared to the data generated bv t he full simulation (for 
further details see lMaris et al.l2005() . A dedicated support 
IDL library to handle the files of expansion coefficients, 
named ZLE_IDL, has been created (jMaris et al...2004b() . 

4. Results 

The contour plot in Fig.Ois an example of Zf{P, R©, Rp) 
for the Smooth component contribution calculated for the 
857 GHz frequency channel assuming Ef = 1. According 
to the discussion in Sect. |21 -E/ « 0.65 and the expected 
observed surface brightnesses are about 2 /3 the values re- 
ported in the plot. The figure represents the variation 
of Zf as a function of the pointing direction P for a 
given combination of R© ad Rp . Contours are drawn for 
Zf = 0.21, 0.22, 0.23, 0.24 0.25, 0.30, . . ., 3.2 MJy/sr. 
Given the cylindrical symmetry of the IDPs, the point- 
ing in the plot is expressed as a function of the ecliptical 
latitude and of the relative ecliptical longitude, i.e. the 
difference between the longitude of the pointing direction 
and the longitude of the solar direction which is at the 
centre of the plot. The blue dotted line represents the 
path described by a Planck beam at the centre of the 



field of view assuming the nominal scanning strategy and 
the grey band represents the region observed by consid- 
ering all the Planck beams. Having the IDPs cloud a 
cylindrical symmetry, and the symmetry reference frame 
being nearly equivalent to the Planck comoving refer- 
ence frame, for the nominal scanning strategy both the 
signal contour levels and the region observed by Planck 
will shift approximately in the same way when the spin 
axis is repointed. Consequently, only a small fraction of 
the possible pointings in the cloud reference frame will 
be observed by Planck. In the case of more complicated 
scanning strategies, such as those including slow spin axis 
precession about the Sun - Satelli te di rection or slow oscil- 
lations above / below the ecliptic l|Dupac fc Taubeii2005j) , 
the scanning path will be shifted normally and along the 
ecliptic plane. Slow spin axis precession or oscillations 
with semi-amplitudes of ~ 5° — 10°, such as those con- 
sidered for the Planck scanning strategy, will change the 
path reported in the plot by ~ 5° — 10° with a result- 
ing signal difference more sensitive near the ecliptic plane 
(~ 0.05 — 0.1 MJy/sr). However for any reasonable scan- 
ning strategy the envelope of all the possible scanning path 
will be only twice or three times wider than the grey band 
in the figure. 

The tilt of the cloud with respect to the ecliptic plane 
and parallactic effects induced by the motion of the space- 
craft with respect to the cloud introduce small modifica- 
tions in the pattern of the contour lines and between the 
signal TODs from different scan circles. 

TODs may be generated at any desired sampling rate, 
for example from that corresponding to 1/3 of the instru- 
mental FWHM resolution to 1° — 2° resolution. It is im- 
portant to estimate the error in the computation of the 
signal in the TODs when the true convolution with the 
beam about its centre direction is replaced by the con- 
volution with a "pencil beam". Denoting with S± and 
(5 1 1 the displacements from the beam centre respectively 
along the direction parallel to the scan circle oriented in 
the scan direction and normal to it towards the Sun, we 
compute the Zf derivatives along these directions. They 
are displayed in Fig. El where dlnZ f /dd± and dlnZf /d5^\ 
are plotted for different positions along the scan circle of 
Fig. Oat 857 GHz and for the Smooth component. For the 
Smooth component \d\nZf/dS±\ < 5.6 x 10"'^ deg"-^ and 
\dlnZf/d6\\ \ < 8.8 x lO'^ deg-^ (here Zf is in MJy/sr). 
For displacements less than FWHM/ 2 ~ 2.5 arcmin the 
error introduced by this approximation is at most 0.03% 
along S± and 0.05% along These error estimates can 
be linearly rescaled to larger displacements S± and ^|| such 
as those associated with the sky pixelization. 

Fig. 01 represents a portion of a TOD, with the asso- 
ciated uncertainty, simulated at 857 GHz without noise 
and with a sampling at a resolution of 1°. It is the gen- 
erated sampling ^smooth. 857 GHz from Fig. O scaled with 
^Smooth, 857 GHz = 0.65. In the figure we report for com- 
parison our preliminary estimate of the sum of the other 
ZLE components. Comparing the TOD with the corre- 
sponding contour map one sees that maxima in ZLE (red 
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G < Gcut A^dsc coy{Z, G) var(Z) var(G) cov(AZ, G) var(AZ) 



[MJy/sr] [MJyVsr^] [MJyVsr^] [MJyVsr^] [MJyVsr^] [MJyVsr^] 



oo 


8010 


1.60E-1 


4.72E-2 


1.12E+3 


3.4 


6.25E-2 


1.87E-3 


6.9E-4 


4.0 


5141 


2.01E-2 


4.62E-2 


9.76E-1 


0.43 


7.57E-3 


1.49E-3 


9.8E-4 


3.0 


4524 


2.38E-2 


4.55E-2 


5.29E-1 


0.52 


4.75E-3 


1.40E-3 


l.lE-3 


2.0 


3591 


1.52E-2 


4.38E-2 


2.12E-1 


0.35 


8.10E-4 


1.18E-3 


1.3E-3 


1.0 


2007 


7.58E-3 


3.27E-2 


3.49E-2 


0.23 


-3.94E-4 


9.09E-4 


2.0E-3 


0.9 


1772 


7.20E-3 


2.95E-2 


2.56E-2 


0.24 


-4.56E-4 


8.67E-4 


2.1E-3 


0.8 


1495 


5.42E-3 


2.48E-2 


1.76E-2 


0.22 


-3.27E-4 


8.09E-4 


2.4E-3 


0.7 


1220 


3.98E-3 


1.95E-2 


1.17E-2 


0.20 


-1.91E-4 


7.68E-4 


2.8E-3 


0.6 


881 


2.32E-3 


1.44E-2 


7.25E-3 


0.16 


-8.37E-5 


6.96E-4 


3.4E-3 


0.5 


505 


1.54E-3 


1.06E-2 


4.17E-3 


0.14 


8.33E-5 


6.67E-4 


4.7E-3 


0.4 


177 


7.58E-4 


5.51E-3 


2.39E-3 


0.14 


9.41E-5 


5.72E-4 


8.5E-3 



Table 2. Correlation between Galaxy, ZLE and Differential ZLE at 857 GHz, for rdsc = 1°, /3 = 85°, (5£ = 2°. 



lines in the figure) occur when the beam crosses the plane 
of the IDPs cloud, slightly below the ecliptic plane for the 
considered case. Minima instead occur when the beam is 
approximately orthogonal to it. For Planck, the spin in- 
duced modulation of the ZLE signal has a main period 
equivalent to 2 cycles per minute (for the nominal spin 
rate of 1 r.p.m.). In a single scan circle, two maxima oc- 
cur when the beam crosses the ascending node and the de- 
scending node between the scan plane and the IDPs cloud 
plane. Since the plane is tilted on the ecliptic. Fig. fallows 
us to predict the occurrence of asymmetries between the 
two peaks even if the scan circle is centred at the antisolar 
direction. The secondary components contribute « 10% 
to the bulk of the ZLE. 

To compare the expected ZLE signal with the Galactic 
foreground, TODs for the Galactic emission have been 
generated and averaged within a circle of 1° radius. 
Galactic emissio n TODs are generat ed using Galactic 
maps ob tained bvlSchlegel et alJ 1 19981) with the prescrip- 
tions in iFinkbeiner et al.l l)l999|) for the scaling of the 
Galactic surface brightness as a function of frequency and 
pointing direction. Due to the tilt of the ZLE symmetry 
plane over the Galactic plane, the ZLE at 857 GHz is com- 
parable to the Galactic emission at low ecliptic latitudes 
where the Galaxy is weak, contributing a peak surface 
brightness of « 0.7^g'2 MJy/sr, or approximately half of 
the weakest Galactic signal along that circle. Of course, 
the ecliptical longitude of the spin axis about which the 
scan circle is drawn will affect the relative contribution of 
the ZLE with respect to the Galaxy. For this reason, we 
report in Fig.|31the variation of the ratio EfZj/Gf (where 
Gf denotes the Galactic surface brightness - we consider 
here only the dominant dust emission) with the ecliptical 
longitude of the spin axis, for three different Planck fre- 
quencies and for circular patches of 1° radius. The white- 
full line represents the ratio averaged over the given scan 
circle. The black band is the ±3tT range of such varia- 
tion, while the highest ratios expected for each scan circle 
are represented by the gray-dashed line. At 857 GHz for 
about half of the scan circles the expected peak ZLE is 
roughly half of the Galactic dust emission. Since the ZLE 
frequency scaling is not much different from that of the 



Galactic dust emission, also at the lower frequencies con- 
sidered here its peak contribution to the sky emission is 
still larger than some ten percent of the Galaxy for most 
of the circles. This contribution is compared to the instru- 
mental sensitivity. The gray band at the bottom of each 
frame in the figure represents the ratio between the in- 
strumental noise, iVy, and the Galactic dust Nf/Gf for 1° 
circular patches, averaged over a scan circle. On average, 
the ZLE contribution is largerly above the instrumental 
noise from 857 GHz to 353 GHz. 

4.1. Time dependence characterization 

The time dependence in the ZLE signal is characterised by 
the short term modulation shown in Fig. 01 and by a long 
term modulation derived from the effective motion of the 
spacecraft within the IDP cloud determined by the L2 or- 
bital motion around the Sun and the spacecraft Lissajous 
orbit around L2. 

The effects of the spacecraft motion are better repre- 
sented in the reference frame defined by the cylindrical 
symmetry of the cloud and corotating with L2 . To under- 
stand the long term modulation we have to consider the 
foUowings: 

1. The tilt of the IDPs fundamental plane over the eclip- 
tic, which introduces a seasonal modulation in the sig- 
nal seen by a spacecraft bound stay near L2; 

2. The ellipticity of the L2 orbit about the Solar System 
barycentre moving the spacecraft with respect to the 
centre of the cloud; 

3. Since the centre of the Smooth component does not 
coincide with the Sun, even a circular orbit around 
the Sun will induce changes in the spacecraft position 
with respect to the cloud centre of symmetry; 

4. Following its Lissajous orbit around L2 the spacecraft 
changes its height over the cloud symmetry plane. 

The spacecraft position is then Rp = Rl2+'5Rlj, with 
Rl2 the position of the L2 point and (5Rlj the Lissajous 
orbit. For typical Planck Lissajous orbits around L2, 
I^RljI is less than a few xlO^ Km. 
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Due to the elhpticity of the Earth orbit, the distance 
of L2 from the Sun varies during the year by a 3%, i.e. of 
0.03 AU ~ 4.5 X 10^ Km. 

Being on the echptic, the L2 point changes its distance 
from the bulk component symmetry plane due to its tilt. 
This induces at maximum a vertical oscillation of ±5.2 x 
10^ Km. In addition the Sun is off-centred with respect 
to the centre of IDP cloud of the Smooth component « 
2 X 10^ km. 

The largest seasonal dependence is due to the tilt of the 
IDP fundamental plane. It affects mainly the value of the 
minima of the surface brightness observed by Planck. 
When the spacecraft is below the fundamental plane of 
the IDP cloud the optical depth towards the North eclip- 
tic Pole is larger than that towards the South ecliptic Pole, 
resulting in a North/South asymmetry in the minima. As 
Planck orbits about the Sun, the spacecraft goes toward 
the node between the ecliptic and the IDP cloud symmetry 
plane, crosses it and enters a region where the symmetry 
plane is below the ecliptic. So, with time the North/South 
asymmetry goes to zero and then reverts its sign. The 
tilt of the symmetry plane over the ecliptic does not sig- 
nificantly affect the level of maximal ZLE surface bright- 
ness observed by Planck, while it affects the location of 
the maxima and the shape about the peak. Fig. |S1 repre- 
sents the modulation of the minima and the North/South 
asymmetry for a 857 GHz horn supposed to be aligned 
with the telescope optical axis. The full curve represents 
the surface brightnesses looking to the North ecliptic Pole. 
The dashed curve represents the surface brightnesses look- 
ing towards the South ecliptic Pole. The relative seasonal 
modulation is about 20%. 

In Fig. m we compare the yearly averaged ZLE surface 
brightness with the daily surface brightness at 857 GHz. 
The variation of the Planck position with respect to 
the symmetry plane of the Smooth component introduces 
variations of up to ~ 10% in the surface brightness with 
respect to the yearly average surface brightness, almost in- 
dependently of the considered frequency and Ef. For typ- 
ical Lissajous orbits, the variation of the Planck height 
with respect to the cloud symmetry plane is ~ 10% of the 
variation induced by the tilt of the symmetry plane on the 
ecliptic but with a periodicity of 6 months and phase dis- 
placement with respect to the yearly periodicity related 
to the exact launch date. Therefore, about 10% of the 
above 10% variations of surface brightness induced by the 
effective Planck motion is introduced by the Lissajous 
orbit. Clearly, this is a second order effect for studies of 
the yearly averaged properties of the Smooth component, 
but it is still larger than the sensitivity of Planck TODs 
averaged over 1° or 2° resolution, as it will be discussed in 
Sect.El Secondary components contribute to about 10% of 
the global ZLE. Therefore, neglecting the Planck orbit 
may significantly reduce the accuracy with which these 
components can be studied. Finally, the differential ap- 
proach to ZLE separation that has several advantages with 
respect to other approaches (see Sect. 15.411 exploits the 
variation of the ZLE during the mission. Neglecting the 




Fig. 7. Absolute and relative variation of the ZLE sur- 
face brightness during the year at 857 GHz. The full-red 
line (light-gray in bw) is the yearly averaged ZLE surface 
brightness [MJy/sr] for a given eclipitical latitude. The 
surrounding green band (dark-gray in bw) is the varia- 
tion of ZLE surface brightness during the year. Below (full 
black lines) the relative variation (with respect to the av- 
erage) of ZLE for a set of representative longitudes (from 
the left of the black curves, from 0° to 75° in steps of 15°, 
respectively from top to bottom) . Note that these relative 
variations reach peaks of ~ 10%. Data are calculated for 
patches of 1° in radius, Es57 ghz = 0.65. 

Lissajous orbit effect will result in an error of 10% in this 
kind of analysis. Of course, the precise inclusion of the 
spacecraft position is not a concern from a computational 
nor practical point of view. 

4.2. Frequency scaling for the spatial distribution 

The frequency scaling for Zf(P,R0,Rp) is a theoretical 
outcome of the model and can be used to check the extent 
by which it is possible to assume Zf(P,R0,Rp) oc f^^''. 
To study the frequency scaling for the spatial distribu- 
tion of ZLE, sets of Zf(P,R0,Rp) have been generated 
for / covering all the Planck frequency channels up to 
/ — 1200 GHz, a fixed combination of Rq and Rp posi- 
tions, Ef = I, and scanning the sky in circles of increasing 
angular radius from 75° to 95° centred in antisolar direc- 
tions. The corresponding data are plotted in the upper 
frame of Fig. |51 while the lower frame represents the spec- 
tral index obtained by fitting a power law dependence for 
the surface brightnesses obtained for a given pointing di- 
rection. 

Note that the average 7/^2 is « 1.971, close to the ex- 
pected value "ff_z — 2. In addition, the spectral index is 
modulated with the pointing and is anticorrelated with the 
surface brightness, the higher spectral indices occurring 
for lower surface brightnesses. The amplitude of the mod- 
ulation, however, is modest, A7/ w 9.4 x 10"'^ in terms of 
peak-to-peak signal. The colour correction would be then 
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CTp 




AEi 


arcmin 


1° Patches 


2° Patches 


0.5 


(6.3 ± 18.4) X 10' 


(1.0 ±2.09) X 10"* 


1.0 


(8.7 ±27.9) X 10" 


(0.9 ± 3.13) X 10"* 


1.5 


(5.2 ±4.10) X 10" 


(1.7 ±4.50) X 10"* 



Table 3. Expected effect at 857 GHz of an isotropic, 
random pointing error on the Ef sensitivity as a function 
of CTp. Calculations are for a cut of 1 MJy/sr on the Galaxy 
and 1° and 2° patches. 



insignificantly affected. The effect of a change of ±0.07 AU 
in the position of Planck over the ecliptic results in a 
^'l'f,z < 0.0027/.2 ^ 0.004 which has negligible effects on 
Kf too. Smaller variations occur shifting the Planck po- 
sition in other directions. The same holds for a shift of the 
Sun position with respect to the centre of the cloud. 

These results assure that a power law is an adequate 
approximation within each frequency band, but this is 
only approximately true considering the full range of 
Planck frequencies. For example, in the range 70 GHz < 
/ < 144 GHz the mean spectral index is 1.990 while in 
the range 545 GHz < f < 1000 GHz the mean spectral 
index is 1.928. 

Assuming at f — 100 GHz the same spectral index as 
at / = 857 GHz, the relative error induced in the sur- 
face brightness prediction at 100 GHz would be « 7%. 
Therefore, a set of spatial templates, one for each fre- 
quency channel, has to be produced in order to simulate 
the spatial dependence of the ZLE with an accuracy of 1%. 
On the contrary, by relaxing the required accuracy to 10%, 
a spatial template at an appropriate reference frequency 
(for example a,t f — 1000 GHz) followed by a spatially 
independent frequency scaling can be used. 

These effects, although very small, may be not neg- 
ligible in CMB studies, that are mainly carried out at 
1/ ~ 100 GHz requiring foreground removal with an accu- 
racy better than ~ 1%. 

5. Separability of ZLE with Planck data 

The ZLE represents a foreground contamination mainly 
relevant for the higher frequency channels of Planck. Its 
contribution has to be removed to accurately study the 
Galactic large scale structure and its frequency scaling. 
Given the weakness of the ZLE signal, the final separation 
quality will may rely on the prior information added to 
the system derived from other missions at IR bands, such 
as IRAS and COBE where the ZLE dominates the sky 
emissivity. We discuss here four different approaches to 
ZLE separation. 

5.1. A test of a "blind" map-based approach 

When ZLE is folded over a map, its histogram is strongly 
non-Gaussian, as for the Galaxy. In principle this sug- 
gests the possibility of obtaining a proper separation of 



the ZLE through blind component separation methods al- 
ready used to analyze microwave maps. 

We performed some numerical experiments with the 
FastICA code (Maino et al. 2002) applied to full-sky maps 
obtained by adding the ZLE and the Galactic emission. 
FastICA is a blind-separation method which usually uses 
as input A'maps maps at different frequencies that are lin- 
ear combinations of iVsignais(< A^maps) signals, all non- 
Gaussian except for at most one. The code gives as out- 
put maps of the various signals. We exploited the two 
frequency channels at 857 GHz and 545 GHz and, as a 
test case, neglected the noise. Although in principle this 
method could be investigated, possibly by applying it to 
Planck channels combined with IR data, this test gives 
discouraging results. Likely this is due to the weakness of 
the ZLE emission and the fact that this code does not use 
any prior information about the ZLE spatial distribution. 
In addition part of this information is lost when passing 
from TODs to maps. 

5.2. A "non-blind" map-based approach 

An "ad hoc" strategy for ZLE detection and separation 
can use a prior information, derived from IR observa- 
tions, on the spatial dependence of the ZLE. We devel- 
oped a model for ZLE detection and separation based 
on the extrapolation of the geometrical information from 
COBE/DIRBE to Planck frequencies (or from any other 
reliable model), leaving as free parameters the emissiv- 
ity corrections, Ef, at Planck frequencies. The ability of 
Planck to measure the ZLE is then translated into its 
accuracy in the determination oi Ef at different frequen- 
cies. 

A non-blind separation based on maps could be dif- 
ferently investigated taking as prior information the ZLE 
spatial dependence and the existing templates of the 
Galactic emission degraded to angular resolutions com- 
parable with the scales of significant ZLE variations (w 
5° - 10°). 

For example, denoting with Sf the map obtained from 
the observed signal minus its average, Qf the Galactic 
emission template (or any other relevant background sig- 
nal) minus its average, Zf the template for ZLE spatial 
distribution minus its average and A// the template for 
noise spatial distribution minus its average for the consid- 
ered scanning strategy and satellite orbit, then the map 
for the signal (minus its average) may be approximated 
with 

Sf^EfZf+afGf+Aff, (15) 

where a/ accounts for possible overall systematic scaling 
errors, both from calibration and frequency extrapolation, 
in the "Galactic" template. It is possible to attempt a 
minimisation of 

= XI - ^f^f.p - "/S/.p) ' (16) 
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where the sum is taken over all the pixel index of the map 
p and cr^ p is the noise variance at each pixel. 

Assuming stationary noise equidistributed over the 
map, the minimisation provides the estimators for Ef 
and af 



A = Tjzz'^gg — (^zg)^ • 



(17) 



Here Ezz = Ep ^.p^/.p/^^L' ^^g = Ep 5/,p6^/,p/f^L' 
SzG = Y.p^f,pQl,p/^),p^ ^zs = Y.pZf,pSf,p/(y)^p, 
^GS = Y.pGf,pSf,p/crjp. 

In addition, it is possible that the true scaling factor 
is not constant over the sky. To simulate this effect we 
replace a constant scaling a j with a normally distributed 
variable, with expectation a/ and RMS Saf. 

We consider here an illustrative case with sensitivity 
per pixel c/^p ~ 4 x 10^'^/ ^JNcSJ^p MJy/sr, where N(,ffj^p 
the effective number of observations made during the mis- 
sion which contributes to the pixel p, Ef = 0.65, a/ = 1.2, 
Ca./ — 0.01 at / = 857 GHz, and a map sampled at ~ 2° 
resolution, and taking all the pixels in the map (i.e. in- 
cluding also regions where the Galaxy largely dominates) 
the RMS for ay is w 2 X 10~^ with a bias of the same 
order, while Ef is recovered with an RMS accuracy of 
about 0.04 but an excess bias of about 0.06. On the other 
hand, removing all the pixels where the signal from the 
Galaxy does not greatly exceed that of ZLE reduces the 
bias. Removing pixels for which the Galaxy exceeds the 
surface brightness of 4 MJy/sr, Ef and af are recovered 
with a RMS accuracy of about 0.01 for af and 0.05 for 
Ef. Their expectations are very close to their input val- 
ues, with biases of a few x 10^'^. (In the remain we consider 
limits of surface brightness at values exceeding 1 MJy/sr). 
As can be seen, the accuracy of the method is very good. 

On the other hand, any (positive or negative) residual 
contribution from the ZLE in the Galactic template from 
the data analysis of the IR data will be scaled to Planck 
frequencies and will introduce a systematic effect which 
will be correlated with the spatial template adopted in 
Eq. (|16|l . This will result in biases in the recovered Ef 
and ctf values. An end-to-end evaluation of this effect, 
beyond the scope of this paper, would require to analyse 
in detail the mission and the data reduction procedure 
used to obtain each IR data set used in preparing the 
Galactic template. However, the results of this approach 
can be compared with the differential method described 
in Sect. l5.4l that automatically by-passes this problem. 

5.3. A "total- power" TODs-based approach 

The prior information derived from IR observations, dis- 
cussed previously also can be applied to the time domain, 
taking also into account the time dependence of the space- 
craft position within the IDP cloud. 

In this approach the separation is based on the knowl- 
edge of the time dependence of the ZLE signal in TODs 



derived from the spatial pattern Zf{P, Hq, Rp). Then, as 
before, we define an estimator Ef of Ef starting from the 
observed data and the known spatial pattern. 

Again, denoting with Gt, Zt, Nt, and St the Galaxy, 
the ZLE, the noise and the signal (minus their average 
values over the mission) , and neglecting systematic instru- 
mental effects 



EfZt + afGt + Nt 



(18) 



This equation is analogous to Eq. H15|l where p is re- 
placed by t and maps are replaced by TODs. So Ef and 
a f can be obtained from Eq. p7|l but replacing p with t 
and maps with TODs. Note that in the case of stationary 
noise uncorrelated with the signal, at is constant all over 
the TOD and substitutions like ^zg — J2t ^t^t/'^t — 
cov{Z,G)/(T'j are allowed ^. The surface in this case 
is similar to the case in Sec. 15.21 but here the sensitivity 
per pixel is constant cr/_t ~ 4 x 10""^ MJy/sr. The final 
sensitivity of this approach is found to be similar to that 
found for the method in Sec. 15.21 

Consider the case in which the Galactic cont r ibutio n 
is neglected in the fitting, as in iKelsall et al.1 l|l998|) . 
Assuming stationary white noise, after minimisation of 



X 



X! 0-2 



EfZt — St 



the estimator formula is 



E 



f 



^zz 



(19) 



(20) 



Assuming that the noise is uncorrelated with the sky 
signal, the expectation for Ef is 



E[Ef] =Ef + 



cov(G,Z) 
var(Z) 



(21) 



This simple estimator is then affected by a bias due to the 
correlation between Z and G. 

The bias is likely negligible for iKelsall et al. I l|l998l) 
since in their case the ZLE signal is much larger than 
the Galaxy emissivity so that |cov(G',Z)| ^ var(Z'). 
But at 857 GHz we obtain cov(G, Z) « 0.16 MJy^/sr^, 
to be compared with var(Z) « 0.047 MJy^/sr^ and 
var(G) w 1117 MJy^sr^, leading to a bias « 3.4 in 
Ef. Such high bias such a small covariance comes from 
the fact that cov(G,Z) is of the same order of magni- 
tude as var(Z). Selection of samples in order to reduce 
cov(G, Z) \ /var(Z) does not mitigate the problem. For ex- 
ample, removing all the samples where the Galaxy is larger 

* We consider in this work instrumental white random noise. 
Planck receivers are affected by l//-like noise that introduces 
long term correlations appearing as offsets in TODs (and as 
stripes in maps). On the other hand, destriping algorithms (see, 
e.g., 'Burigana ct al. 1997) accurately remove this effect in the 
TODs (and in the maps). 

^ This is not true for a map since in this case is a function 
of p even for stationary white noise. 
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Freq 


Minimum AEf 


Optimal Gcut 


Min. Gcut 


Max. Gcut 


[GHz] 




[MJy/sr] 


[MJy/sr] 


[MJy/sr] 


857 


0.022 


6.5 


2.3 


26 


545 


0.017 


2.0 


1.0 


8.3 


353 


0.011 


0.50 


0.20 


2.0 



Table 4. Effect of a random relative calibration error E,-^2° = 1% on the Ef determination at different Planck 
frequencies. Column 2 is the minimum RMS error for the optimal surface brightness cut reported in Column 3. 
Columns 4 and 5 give the range for which the RMS error is less than twice the value in Column 2. 



than 1 MJy/sr leads to cov(G', Z) « 7.6 x lO^^ MJyVsr^ 
var(Z) w 0.033 MJyVsr^, var(G') w 0.034 MJyVsr^ 
with a bias w 0.23. Removing all samples where the 
Galaxy is larger than 0.4 MJy/sr leads to cov(G,Z) « 
7.6 X 10-" MJyVsr^, var(Z) « 5.5 x lO'^ MJyVsr^ 
var(G) « 2.4 x lO'^ MJy^/sr^, with a bias « 0.14. As 
evident, the bias decreases when applying stronger cuts 
but it still remains significant (see Table |2 columns 1 to 
6). 

5.4. A "differential" TO Ds- based approach 

Planck will scan the sky at least twice during the mis- 
sion. Therefore, most of the sky directions will be observed 
at least twice with Planck in different positions within 
the IDP cloud. In the ideal case, the difference between 
these two measures will be due to the difference in the 
ZLE contribution that can be predicted from our model 
plus noise. 

We denote with ti{P) and iii(P) the epochs of the 
first and the second observation of a region seen in the 
direction P and with ^i, Su {Zi, Z\\ or Gi, G\\ or iVi, N\\) 
the corresponding observed surface brightness total signal 
variations (ZLE spatial distribution or Galactic emission 
or noise) with respect to the mean. With these definitions 
the differential surface brightnesses will be 



A^t = ^tjj - St^ = AGt + S/AZt + ATVt 



(22) 



with the convention that i = t\ and that AG* — G\\ — G\ 
and so on. The is defined now as 



Ej^Zt ~ ASt 



(23) 



'AN,t 



where a^N.t is the RMS of the noise for the considered 
samples (for stationary noise aAN,t = "'^^ + '''ajv tu ~ 
an) giving 



2a2 



E 



f 



AZ 



^AZ.AZ 



(24) 



where T^as.az = J2t ^^tAZt/aj^j^ ^ oc cov(A5, AZ) and 
Eaz,az — /'^AN t ^ var(AZ). The expectation 

for this estimator is 



E[Ef] = Ef + 



cov(AG, AZ) _ 
var(AZ) ' 



(25) 



since, by definition, P(ti) 
any pair of ti, tu giving 

E[Ef]^Ef, 



P(iii), we have AG = for 



(26) 



without any bias and without the need to use Galactic 
templates. 

Applying the standard error propagation formula to 
Eq. (|26|l and considering Eq. 124() . if the noise ANt can 
be approximated as stationary, independent and Gaussian 
with variance (tan\ = '^an'^, summing over all of the 
pairs we have 



var 



^AS,AZ 
^AZ,AZ 



j:,n2ANt^]AZt 



J2t '^AN.t 



AZt^ 



Sa 



Z,AZ 



Sa 



Z,AZ 



After some algebra we have 



va,T[E] 



2a an'^ 



7Vcpivar(AZ) ' 



(27) 



where iVcpi is the number of independent (ii,iii) pairs. 
Note that var[i5] does not significantly depend on the 
adopted radius of the patch because of the invariance of 
the CTAAf^Z-ZVcpi ratio. 

Fig. ^2 represents the expected EfAZ for Ef — 0.65, 
the nominal scanning strategy of Planck at 857 GHz, and 
a set of selected ecliptical longitudes of Planck ^. Note 
that the peak differential signal is about 10° — 20° from 
the ecliptic plane. The signal is calculated averaging over 
independent circular patches of 2° in radius. The differ- 
ence between the first and the second scan never exceeds 
0.06 MJy/sr, then being a » 10% effect. However, when 
compared to the sensitivity expected in this frequency 
channel (represented at ba by the gray band), this signal is 
clearly detectable, particularly when the spacecraft is lo- 
cated at ecliptical longitudes of ~ 45° and ~ 270°, where a 
particularly good peak S/N{c± 40) is expected. For a S/N 
threshold « 5, a clear detection of the differential signal is 
expected for Ef > 0.06. In case oiEf « 0.2-0.5, the S/N 
ratio is so good as to open the possibihty to also improve 
the parameters of the geometrical model, to study possi- 
ble spatial dependences in E'/, and to identify secondary 
components. 

® In this calculation it is assumed that Planck is orbiting 
about the L2 point according to the nominal orbit expected 
for a launch in February 2007. Changes in this orbit will only 
slightly change the results discussed above. 
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The noise statistics "per patch" (RMS and patch-to- 
patch correlation) depend on the method used to assemble 
samples from TODs to form patches of sky. To determine 
a simple statistic for the noise, we construct patches i) 
of fixed solid angle (namely, circular patches with radius, 
Tdsc, of 1° or 2°), ii) observed in both surveys with a sig- 
nificantly large number of samples coadded so as to avoid 
a significant difference in the effective weight of each sky 
direction in the two surveys (in order to assure a simi- 
lar coverage of the same patch in the two surveys and 
smooth out possible particularly bright pixels - see also 
the discussion in Sect. 15.5.1(1 . iii) taken contiguously in 
time {At < 1 day x (rdsc/l")) in each survey, and iv) 
avoiding the presence of overlapping patches. These con- 
straints only slightly reduce the number of samples used 
in the analysis. After two surveys the average instrumen- 
tal noise RMS on a single squared pixel with side equal 
to 6 = 5 arcmin for the reference frequency channel at 
857 GHz is 43 mJy. Composing these pixels to form cir- 
cular patches of radius rdsc, the noise per patch is 

fjAAT ~ 43 X lO-^MJy X = 

^ MJy/sr , 

where the factor 2 comes from the fact that we are con- 
sidering the difference between the values observed in the 
same patch in the two surveys taken separately (the vari- 
ance of samples entering in the patch can be neglected 
according to condition ii)). 

Columns 7 to 9 of Table [5] shows the statistics for 
the simulated scan at 857 GHz. The last column gives 
the expected 1 a error on Ef determination according to 
Eq. ((77|l with the expected level of noise. 

Accepting all the sky samples at 857 GHz a RMS ac- 
curacy aEf ~ 7 X 10"'' should be expected. The effect of 
cuts based on the Galactic surface brightness is shown in 
Table 121 By considering regions where the Galactic signal 
is smaller than Gcut = 1 MJy/sr (Gcut = 0.4 MJy/sr) the 
accuracy reduces to aEf ~ 2 x 10~^ {asf 8 x 10"'^). 
On the other hand, for Gcut = 4 MJy/sr (relevant to re- 
duce the impact of the relative calibration uncertainty, 
see Sect 15. 5."^ the number of independent pairs decreases 
from ~ 8 x 10^ to ~ 5.1 x 10^ with ^ 10"^, only 
slightly degraded with respect to the full sky analysis. 
The cut changes the sign of column 7 due to the tilt of the 
ecliptic plane, so that the ZLE and its variation is stronger 
where the Galaxy is weaker. 

A similar analysis carried out at 545 GHz and 353 GHz 
gives analogous results on the scaling introduced by the 
surface brightness cut. By considering the sky regions 
identified by the Gcut = 4 MJy/sr at 857 GHz, with 
the same kind of analysis we find ^ 2.1 x 10^'' and 
aEf ^ 2.6 X 10~3 at 545 GHz and 353 GHz, respectively. 

5.5. Systematic effects in the differential approach 

The most important source of error in the determination 
of i?/ is the bias induced by cov(AG, AZ) when it is com- 



parable to var(AZ). In this subsection the main sources 
of this bias are discussed. 

5.5.1. Error Induced by Sky Sampling 

In the differential approach presented in the previous sec- 
tion we cancel out the Galactic signal. This is a good ap- 
proximation provided that the patch is equally sampled 
in each of the two surveys, implying AG(P) = 0. In real- 
ity, uncertainties in the spacecraft attitude reconstruction 
and in pointing maneuvers and, in particular, the displace- 
ments between the positions of the various samples taken 
in the two surveys will imply that the same patch is sam- 
pled in a different manner in each survey. An estimate 
of the maximum displacement between the position of a 
given sample in the first survey and the nearest sample 
taken in the second survey is given by the maximum be- 
tween half of the spin axis displacement (1.25 arcmin) 
and half of the angular sampling along the scan circle 
(- 0.5 FWHM/3, equivalent to - 0.8 arcmin at 857 GHz). 
The pointing accuracy for Planck is expected t o be bet- 
ter t han .5 a rcmin (la) for each sample ( Burig ana et alJ 
l200ltlPuget et alJ200l . a value smaller than the above es- 
timate ^. Denoting with P the pointing direction neglect- 
ing the displacements discussed above and with Pi(P) the 
effective pointing direction in the i-th survey, the displace- 
ment is (SPi(P) = P,(P)-P. We are interested in AG(P) 
which depends on the combined displacements AP(P) = 
5Pii(P) - <5Pi(P). Then AG(P) ~ VG(P) • AP(P) / 0, 
leading to a bias in the determination of Ef of the order 
of cov(VG • AP^AZ)/var(AZ). 

If we assume that 5Pi does not correlate with the 
Galactic emission or the ZLE, as is reasonable, no bias 
will be introduced. On the other hand, this effect will de- 
grade the sensitivity by 

D,G,patch,2 fOQ'\ 
"Ef. point = (^^(AZ,)2)2 ' ^'''> 

where the variance cr^ q p^^^^j^ is taken over the consid- 
ered patch. We assume that SPi is a bivariate random 
variable, normally distributed, with null expectation and 
isotropic covariance matrix Cp = (Jp~^diag(l, 1), with the 
displacement variance CTp^ 1 arcmin^ uniform all over 
the sky. Therefore, in the above linear approximation, for 
each patch of iVsmp samples the averaged surface bright- 
ness difference will be normally distributed with null ex- 
pectation and variance 

2 

0-D,G,patch ~ varpatch(VG)crp^ , (29) 

*smp 

where the variance is taken over the considered patch and 
the factor 2 is due to the differential approach. For patches 
of 1° (or 2°) radius, and samples of 5/3 arcmin along the 

^ In addition, while finalizing this paper, a significant im- 
provement in the Planck star trackers was achieved lead- 
ing to an expected po inting accuracy of few arcsec (see, e.g., 
iHarrison &: van Leeuwen. 1.2005.'l 'l. 
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Fig. 8. Seasonal modulation of the Smooth component of 
the ZLE at 857 GHz for four selected pointing directions: 
the North ecliptic Pole (N), the South eclitptic Pole (S), 
the Forward direction with respect to the Planck motion 
(FW) and the Backward direction (BW). Longitudes are 
relative to the longitude of the ascending node. Upper 
frame: the surface brightness variation looking to N (full 
line), S (long dashed line), FW (short dashed line), BW 
(dot-dashed line). The variation is the difference between 
the surface brightness along the direction of choice and 
its yearly average. Here i^Smooth, 857 GHz = 0.65. Lower 
frame: North - South (full line) and Forward - Backward 
(dashed line) asymmetries. Asymmetries are defined as 
A = 2{F2 - Fi)/{F2 + Fi) with F2 = ^n, Fi = Fg for the 
North - South asymmetry and F2 = i^pw, ^1 ~ Fbw in 
the other case. 

scan circle and 5/3 arcmin transversally to it, A^smp ~ 
2.5 X 10^ (or « 1 X lO'*). For a simple determination we 
can assume that varpatch(VG') w varpatch(G)/,^sjjjp where 
i?sjjjp is the typical solid angle of the sample, so that 

2 E»varpatch,»(G)(AZO^ a/ 2 

''Ef. ^>oint- (E.(A^.)2)2 e'smpN^mp ■ ^ ^ 

From the simulations of Sect. 15.41 for patches of 1° (or 
2°), cTp = 1 arcmin and applying a cut to remove patches 
where the Galaxy is stronger than 1 MJy/sr, we estimated 
a aj^^ p^jj^j of 6 X 10^'* (or 2 x 10^''). Wc simulate in de- 
tail the effect of a random pointing modelled as above. 
We obtain that, apart from a few spikes that are eas- 
ily filtered, the perturbation induced by the Galaxy is in 
general small. The total error is dominated by the bias 
until strong cuts are applied. The total error does not 
change with Ef. The total pointing error is small. The re- 
sults are shown in Tabled reporting both the expectation 
and the variance of Ai?j. The expectation is small while 
the RMS of the Ef error is consistent with Eq. The 
RMS from our MonteCarlo simulation scales quite well 
with y'2/A^s„ipO'p/£smp- In summary, the random point- 
ing error does not seriously limit the recovery of Ef. 



We consider here the requirement ii) of Sect. 15.41 
We could have a different number of samples ANgmp 
(<?C Ngynp) in the same patch in the two surveys both 
because of the different samples at the boundary of the 
patch and because of the result of the effective scanning 
strategy. This implies a difference in the average Galactic 
signals obtained in the two surveys, related to the fluc- 
tuations of the Galaxy within the patch. With a simple 
algebra it is straightforward to derive that in this case the 
variance of the induced AG in our differential approach is 

^B.G.patch varpatch(Cx)— . (61) 

J^smp 

Eqs. (|29|l and H31() clearly show the relevance of the con- 
straint ii) of Sect. 15.41 in the construction of the patches. 
An approximate comparison between Eq. H29|l and 

Eq. (EU gives G.patcJ'^S.G.patch ~ CpAiVsmp/2ap2. 

As 4mp > cFp, we have al^^p^^^,^ > cr^ ^^^^^^^ even 
for AA'smp ~ 1. This implies that the effect associated 
with AA'smp could significantly impact the final result, if 
not taken properly into account in the data analysis. Of 
course, it has no physical meaning to have a sensitivity 
degradation in the presence of an increased number of ob- 
servations in one of the two surveys. In reality, the above 
computation underlines the relevance of properly assem- 
bling the "elementary" samples in the two surveys in order 
to count the signal from the same sky direction the same 
number of times in each of the two averages over the con- 
sidered patch. 

The worst effect of systematic pointing errors would 
occur in sky regions with bright point-sources. It would 
then be preferable to remove pixels affected by bright 
sources before of the computation of the averages of the 
signals in each patch in order to manage only with signals 
dominated by the diffuse components. 

An example of an intrinsic source of systematic point- 
ing errors is the aberration of light. When not ac- 
counted for, the aberration due to satellite motion about 
the Sun may induce at most a pointing error of « 
2?;/c rad ~ 0.7 arcmin, dominated by the Earth mo- 
tion about the Sun. The factor of 2 comes from the fact 
that patches are acquired at most at about 180° of longi- 
tude when the orbital motions are toward opposed direc- 
tions in the sky. This effect may introduce an important 
bias; however since the spacecraft velocity is known within 
1 Km/sec or better, the effect may be removed by correct- 
ing the selected pointings. 

5.5.2. Doppler shift 

The relative motion of the satellite with respect to the 
Sun induces a Doppler shift in the Galactic signal ob- 
served during each of the two scans. The effect will be 
^f/f ~ 2v/c X lO"** and assuming G oc /" the surface 
brightness variation will be 6G/G ~ a2v/c. Even for a in 
the range 2-3, \SG/G\ < 8 x 10"'*. Assuming AG w 5G 
and taking the statistics from the signal correlations from 
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Fig. 9. Deviation of the ZLE from an scahng as a func- 
tion of the direction in the sky. Upper frame: the original 
data for / = 100, 143, 217, 353, 545, 857 GHz (from lower 
to upper surface brightness). Lower frame: deviation from 
the scaling in terms of 2 — 7/,z) where 7/^^ is obtained 
by fitting a power law scaling on the data in the upper 
frame. In the horizontal axis we report the pixel indices 
for an HEALPix map with A^sido = 16 ordered according to 
the ring scheme. The colatitude in the map decreases from 
left to right. (Colour figure on the electronic version.) 



TablcH then cov(AG, AZ) < a2(t)/c)cov(G, AZ) w 1.6 x 
IQ-s MJy2/sr2 equivalent to a bias jAS/l < 2.6 x 10"^^ 
However the application of 4 MJy/sr cut will reduce this 
bias by an order of magnitude. Of course, a further rele- 
vant reduction (by a factor ~ 10 or 100) of this effect can 
be reached with a simple modelling of the Galatic emission 
spectral index (for example at ^ 10% or 1% accuracy). 

In addition, assuming the IDP cloud to be at rest 
around the Sun (neither shifting, nor rotating in time), 
the Dopplcr shift will affect the ZLE surface brightness 
too, so that a \EfAZ\ < EfZ2v/c « 1.8 x IQ-* MJy/sr 

equivalent to a bias in Ef of at most 10^*. 
5.5.3. Calibration uncertainty 

The impact of any absolute calibration error of Planck 
data uniform all over the mission is of simple evaluation. 
In the differential approach, any calibration uncertainty 
of this kind produces a null effect in the Galactic signal. 
The only final error will be a multiplicative uncertainty in 
the ZLE and then in Ef given by Ai?/^g — Ef \dg\, where 
5g is the relative uncertainty of the absolute calibration 
(1% — 3% according to FIRAS absolute calibration accu- 
racy from - 300 GHz to - 900 GHz). 

The evaluation of the impact of calibration errors vary- 
ing during the mission is more difficult (i.e. relative cal- 
ibration errors). Clearly, the difference in the calibration 
errors in each patch in the two considered surveys pro- 



Fig. 10. Map of ^Smooth, 857 GHz with A^sidc = 128. The 

top of the map corresponds to the direction of the anticen- 
ter of the cloud. The sides and the middle of the map to 
poiting directions in the cloud symmetry plane. Compared 
to Fig. 121 the circles phase is 0° when the pointing direc- 
tion lies on the equator near Xcomt and 90° on the equator 
near Fcorot (the North pole of the cloud). (Colour figure 
on the electronic version.) 



duces systematic effects in our differential approach pro- 
portional to 6g cov(G, AZ)/var(AZ). 

We evaluated the implications of random relative 
calibration errors. The relative calibration accuracy of 
Planck/HFI highest frequency c hannels is currently un- 
der definition. IPaiot et"al] l(2000l) (§ 3.3.6.1) reported a 
prelim inary relat i ve pix el-to-pixel calibration accuracy of 
3%. |Piat et"al] l)2002|) investigated HFI calibration with 
a kinematic dipole at frequencies at ly ^ 400 GHz and 
with Galactic templates at ^ 400 GHz. The authors re- 
ported a relative calibration accuracy on each scan circle 
of ~ 10% at 545 GHz and of - 3% at 353 GHz with 
suitable choices of the surface brightn ess cut. 

Rescaling the error estimate by IPaiot et al.l l(2000|) 
from ~ 5 arcmin pixels to circular patches of ~ 2° radius, 
w e find a relative a ccuracy, i?r,2°, of 0.06%. The results 
of IPiat et all ll2002l) can be rescaled to the above circular 
patches considering that we have about 10^ scan circles 
per circular patch. In this way we estimate i?r,2° ~ 1% at 
545 GHz and Er,2'' 0.3% at 353 GHz. 

We then consider i?r,2° in the range 0.1% — 1%, so per- 
turbing the simulated signal from each patch. The com- 
parison between the estimation oi Ef obtained in the ab- 
sence of this systematic error and by repeating the per- 
turbed simulation described above for many thousands 
of realizations and for various Galactic cuts is shown in 
Fig. El f^or a particular case. As evident, at 857 GHz, 
assuming i?r,2° — 1% and Galactic cuts in the range 
Gcut '-2-26 MJy/sr we find a RMS error on Ef of 
^ 3 X 10^^ (with an expectation value of the residual bias 
of < 10-4). 
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Fig. 11. Differential surface brightness of the ZLE from 
the Smooth component in the 857 GHz channel, Ef — 0.65 
calculated for a spacecraft ecliptical longitudes 77.7° (solid 
line), 122.7° (dotted line), 167.7° (dashed line), 212.7° 
(dash dot line), 257.7° (long dashed hue), 302.7° (solid 
thick line), 347.7° (dotted thick line), 392.7° (dashed thick 
line). The ascending node of the symmetry plane for the 
Smooth component cloud is assumed to be at a longitude 
of 77.7°, so that these longitudes are equivalent to an an- 
gular distance of of 0°, 45°, 90°, 135°, 180°, 225°, 270°, 
315° from the node respectively. The surface brightness 
scale is in units of MJy/sr. The yellow (light-gray in bw) 
band represents the 5a sensitivity (white noise limited) 
expected from Planck in the 857 GHz frequency channel 
for circular patches of 2° of radius. (Colour figure on the 
electronic version.) 

For £'r,2° ^ 1% the RMS errors on Ef are - 2.5 x 10"^ 
or 1.5 X 10^^ respectively at 545 GHz or at 353 GHz 
for suitable Galactic cuts. Table ^summarizes our results. 
The range of surface brightness cuts appropriate for each 
frequency allowing the above sensitivity spans about one 
order of magnitude. We verify that for £'r.2° 5% the 
RMS error on Ef scales linearly with £'r,2° at all frequen- 
cies. 

Of course, the optimal Galactic surface brightness cut 
should be identified a posteriori according to the recov- 
ered value oi Ef and with a proper trade-off between the 
various systematic errors. 



6. Conclusions 

We presented an analysis to predict the level of contam- 
ination from the Zodiacal Light Emission (ZLE) in the 
survey of the forthcoming Planck mission and to evalu- 
ate the ability to extract the ZLE signal from the Planck 
data. This signal is used to gain more information about 
the ZLE physical properties. 

Our starting point is the model of lKelsall et al. I 



1 

for the ZLE based on the COBE data describing in de- 
tail the emissivity of the IDP cloud for wavelengths up 



to about 300 um con iplemented by the recent result of 
iFixsen fc Dwekl l|2002f l on the spectral behaviour of the 
ZLE. According to the COBE model, four components 
contribute to the ZLE, the dominating smooth compo- 
nents on which this paper is focussed, the Earth orbit 
locked ring of dust (or circumsolar ring), the trailing blob, 
and three bands of dust. 

With respect to other foregrounds usually considered 
in CMB studies, the ZLE (as the other Solar System ob- 
jects) is peculiar, its surface brightness depending not only 
on the pointing direction but also on the instantaneous 
position of the observer within the Solar System. This un- 
derlines the relevance of a study of the ZLE not only on 
maps but also on time ordered data streams. Taking the 
average yearly position of Planck in the IDP dust cloud 
will result in a ~ 10% error in the ZLE surface bright- 
ness estimate, while neglecting the Planck orbit about 
L2 will imply a ^ 1% error in the ZLE surface brightness 
estimate. Since the ZLE differential approach separation 
exploits the ~ 10% variation of the ZLE surface brightness 
between the two surveys, neglecting the Planck orbit will 
introduce a non-negligible systematic ~ 10% error in this 
method. This will have a large impact on the study of 
secondary ZLE components, a natural extension of this 
work. 

We have implemented the COBE model in a dedicated 
program that computes Zf(P, Rq, Rp) for a given list of 
values of P, R©, Rp and a set of parameters describ- 
ing the properties of the ZLE component for which the 
calculation has to be performed. Of them, the emissivity 
correction Ef for each component and for each frequency 
band is the hardest to extrapolate to Planck frequencies 
and it carries most of the physical information on the IDP 
population producing the bulk emission at these frequency 
bands. We then focused on the capability of Planck to 
recover Ef. Since the Ef scalings are largely uncertain, the 
code separately generates the surface brightnesses appro- 
priate for each desired component. The outputs produced 
for the various components then can be easily combined 
according to the user need. 

Since the TODs are large and their realization is 
strictly related to the effective scanning strategy and satel- 
lite orbit we also implemented an approximate, but accu- 
rate, method able to compress the R0, Rp dependencies 
in TODs for the desired pointings in matrices of appropri- 
ate series expansion coefficients. 

We investigated the possibility of obtaining a proper 
separation of the ZLE through a blind component separa- 
tion method (FastICA) already used to analyse microwave 
maps, without obtaining encouraging results. 

A non-blind separation based on maps has been inves- 
tigated taking as prior information the ZLE spatial de- 
pendence and the existing templates of the Galactic emis- 
sion degraded to angular resolutions comparable with the 
scales of significant ZLE variations. We considered here in 
detail the case of the HFI 857 GHz channel. In this case 
the main foreground is the dust Galactic emission. For ap- 
propriate Galactic surface brightness cuts 1 MJy/sr) 
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the value oi Ef recovered with this approach is in good 
agreement with the input one with an RMS absolute ac- 
curacy ~ 0.05. On the other hand, any possible residual 
contribution from the ZLE left in the IR Galactic tem- 
plate adopted from their original data analysis procedure 
will be scaled to Planck frequencies and will introduce 
a systematic effect which will be correlated with the spa- 
tial template adopted in ZLE estimates, possibly result- 
ing in biases in the recovered Ef. A total power approach 
on TOD has been also investigated, providing similar re- 
sults. These two methods require the use of low resolution 
Galactic templates. 

To circumvent the need for Galactic templates and to 
better take into account the effects introduced by the vary- 
ing spacecraft position, we have studied a differential ap- 
proach exploiting the fact that Planck will scan the sky 
at least twice during the mission, staying in different posi- 
tions within the IDP cloud. In the ideal case, the difference 
between these two measures will be due to the difference 
in the ZLE contribution that can be predicted from our 
model plus noise. We find a typical absolute RMS uncer- 
tainty on Ef induced by the limited instrumental sensitiv- 
ity of ^ 10~^ (Icr) at 857 GHz for suitable choices of sur- 
face brightness cuts (or 2.1 x ~ 10"'^ and 2.6x ~ 10~^ at 
545 GHz and 353 GHz, respectively). For typical expected 
values of Ef (« 0.65, 0.26, 0.11 for 857 GHz, 545 GHz, 
353 GHz) the Planck sensitivity will allow an Ef recov- 
ery at 0.15%, 0.8% and 2.4% (Icr) accuracy at 857 GHz, 
545 GHz and 353 GHz, respectively. 

We investigated the impact of the most relevant sys- 
tematic effects, pointing and sampling uncertainty, aber- 
ration of light, Doppler shift and relative calibration un- 
certainty. While the first four effects are expected to be 
significantly below the noise, the last effect is potentially 
the most critical one. For a relative calibration error RMS 
of ~ 1% (0.1%) on patches of 2° radius, we find an abso- 
lute RMS error on Ef of - 0.01 - 0.04 (- 0.001 - 0.004) 
with only a weak dependence on the frequency in the 
range 300 — 900 GHz, corresponding to relative errors 
on Ef - 4%, 10%, 23% (- 0.4%, 1%, 2%) respectively 
at 857 GHz, 545 GHz, 353 GHz for the most hkely Ef 
values expected on the basis of COBE/FIRAS data. This 
may ultimately determine the final accuracy of the ZLE 
extraction from Planck data. 

A web page containing documentation, FITS files, and 
IDL routines to handle the series expansion for ZLE sim- 
ulations is in preparation. 
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